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1. INTRODUCTION 

Both granular superconductors and artificially prepared arrays of Josephson junc- 
tions (JJAs) proved useful in studying the numerous quantum (charging) effects in 
these interesting systems, including blockade of Cooper pair tunneling |I] , Bloch oscil- 
lations |2j, propagation of quantum ballistic vortices jSj, spin-tunneling related effects 
using specially designed S'FS'-type junctions [HE], novel Coulomb effects in SINIS- 
type nanoscale junctions dynamical AC reentrance [Zj and geometric quantization 
phenomena jH] in 2D JJAs (see, e.g., Ref. [Oj for the recent review on charge and spin 
effects in mesoscopic 2D Josephson junctions). 

More recently, it was realized that JJAs can be also used as quantum channels to 
transfer quantum information between distant sites fUJ El E] through the implemen- 
tation of the so-called superconducting qubits which take advantage of both charge and 
phase degrees of freedom (see, e.g., Ref. [13] for a review on quantum-state engineering 
with Josephson-j unction devices). 

Recent imaging of the granular structure in underdoped Bi 2 Sr 2 CaCu20s+s crys- 
tals [Hj, revealed an apparent segregation of its electronic structure into superconduct- 
ing domains (of the order of a few nanometers) located in an electronically distinct 
background. In particular, it was found that at low levels of hole doping (5 < 0.2), 
the holes become concentrated at certain hole-rich domains. Tunneling between such 
domains leads to intrinsic granular superconductivity (GS) in high-T c superconductors 
(HTS). Probably one of the first examples of GS was observed in Y BaiCu^O-j-?, sin- 
gle crystals in the form of the so-called "fishtail" anomaly of magnetization 15 . The 
granular behavior has been related to the 2D clusters of oxygen defects forming twin 
boundaries (TBs) or dislocation walls within CuO plane that restrict supercurrent flow 
and allow excess flux to enter the crystal. Indeed, there are serious arguments to con- 
sider the TB in HTS as insulating regions of the Josephson SIS-type structure. An 
average distance between boundaries is essentially less than the grain size. In particu- 
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lar, the networks of localized grain boundary dislocations with the spacing ranged from 
lOnm to lOOnm have been observed ^3] which produce effectively continuous normal 
or insulating barriers at the grain boundaries. It was also verified that the processes of 
the oxygen ordering in HTS leads to the continuous change of the lattice period along 
TB with the change of the oxygen content. Besides, a destruction of bulk supercon- 
ductivity in these non-stoichiometric materials with increasing the oxygen deficiency 
parameter 5 was found to follow a classical percolation theory [TB] . 

In addition to their importance for understanding the underlying microscopic mech- 
anisms governing HTS materials, the above experiments can provide rather versatile 
tools for designing chemically-controlled atomic scale Josephson junctions (JJs) and 
their arrays (JJAs) with pre-selected properties needed for manufacturing the modern 
quantum devices [TUUTT]. Moreover, as we shall see below, GS based phenomena can 
shed some light on the origin and evolution of the so-called paramagnetic Meissner effect 
(PME) which manifests itself both in high-T c and conventional superconductors [T^ITTI] 
and is usually associated with the presence of 7r-junctions and/or unconventional (d- 
wave) pairing symmetry. 

This Chapter reviews some of the recently suggested novel effects which are ex- 
pected to occur in intrinsically granular non-stoichiometric material modeled by 2D 
JJAs which are created by a regular 2D network of twin-boundary (TB) disloca- 
tions with strain fields acting as an insulating barrier between hole-rich domains in 
underdoped crystals. In Section 2 we consider phase-related magnetization effects, 
including Josephson chemomagnetism (chemically induced magnetic moment in zero 
applied magnetic field) and its influence on a low-field magnetization (chemically in- 
duced PME), and magnetoconcentration effect (creation of extra oxygen vacancies in 
applied magnetic field) and its influence on a high-field magnetization (chemically in- 
duced analog of "fishtail" anomaly). Section 3 addresses charge-related phenomena 
which are actually dual to the chemomagnetic effects described in Section 2. More 
specifically, we discuss a possible existence of a non-zero electric polarization (chemo- 
magnetoelectic effect) and the related change of the charge balance in intrinsically 
granular non-stoichiometric material under the influence of an applied magnetic field. 
In particular, we predict an anomalous low-field magnetic behavior of the effective 
junction charge and concomitant magnetocapacitance in paramagnetic Meissner phase 
and a charge analog of "fishtail" -like anomaly at high magnetic fields as well as field- 
dependent weakening of the chemically-induced Coulomb blockade. 

2. CHEMOMAGNETISM AND FISHTAIL ANOMALY 

As is well-known, the presence of a homogeneous chemical potential fi through 
a single JJ leads to the AC Josephson effect with time dependent phase difference 
d<fi/dt = fi/h. In this paper, we will consider some effects in dislocation induced 
JJ caused by a local variation of excess hole concentration c(x) under the chemical 
pressure (described by inhomogeneous chemical potential /i(x)) equivalent to presence 
of the strain field of 2D dislocation array e(x) forming this Josephson contact. 
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To understand how GS manifests itself in non-stoichiometric crystals, let us in- 
voke an analogy with the previously discussed dislocation models of twinning-induced 
superconductivity [2H] and grain-boundary Josephson junctions [21]. Recall that un- 
der plastic deformation, grain boundaries (GBs) (which are the natural sources of weak 
links in HTS), move rather rapidly via the movement of the grain boundary dislocations 
(GBDs) comprising these GBs. Using the above evidence, in the previous paper [21] we 
studied numerous piezomagnetic effects in granular superconductors under mechanical 
loading. At the same time, observed [HI HH1 [221 [231 Ell m HTS single crystals regular 
2D dislocation networks of oxygen depleted regions (generated by the dissociation of 
< 110 > twinning dislocations) with the size do of a few Burgers vectors, forming a 
triangular lattice with a spacing d > do ranging from lOnm to lOOnm, can provide 
quite a realistic possibility for existence of 2D Josephson network within CuO plane. 
Recall furthermore that in a c?-wave orthorhombic YBCO crystal TBs are represented 
by tetragonal regions (in which all dislocations are equally spaced by do and have the 
same Burgers vector a parallel to j/-axis within CuO plane) which produce screened 
strain fields [23] e(x) = e e~' x '/ d ° with | x | = \/x 2 + y 2 . 

Though in Y Bci2Cuj,Ot_$ the ordinary oxygen diffusion D = Doe~ Ud l kBT is ex- 
tremely slow even near T c (due to a rather high value of the activation energy Ud in 
these materials, typically Ud — leV), in underdoped crystals (with oxygen-induced 
dislocations) there is a real possibility to facilitate oxygen transport via the so-called 
osmotic (pumping) mechanism [23 EE] which relates a local value of the chemical po- 
tential (chemical pressure) /x(x) = /x(0) + V/x ■ x with a local concentration of point 
defects as follows c(x) = e~ M( - x ^ fcsT . Indeed, when in such a crystal there exists a 
nonequilibrium concentration of vacancies, dislocation is moved for atomic distance 
a by adding excess vacancies to the extraplane edge. The produced work is simply 
equal to the chemical potential of added vacancies. What is important, this mech- 
anism allows us to explicitly incorporate the oxygen deficiency parameter 5 into our 
model by relating it to the excess oxygen concentration of vacancies c v = c(0) as 
follows 5 = 1 — c v . As a result, the chemical potential of the single vacancy reads 
u v = /i(0) = —ksT log(l — 5) ~ ksTS. Remarkably, the same osmotic mechanism 
was used by Gurevich and Pashitskii [23] to discuss the modification of oxygen va- 
cancies concentration in the presence of the TB strain field. In particular, they argue 
that the change of e(x) under an applied or chemically induced pressure results in a 
significant oxygen redistribution producing a highly inhomogeneous filamentary struc- 
ture of oxygen-deficient nonsuperconducting regions along GB P4 3 (for underdoped 
superconductors, the vacancies tend to concentrate in the regions of compressed ma- 
terial). Hence, assuming the following connection between the variation of mechanical 
and chemical properties of planar defects, namely /x(x) = KQ e {' x ) (where Qq is an 
effective atomic volume of the vacancy and K is the bulk elastic modulus), we can 
study the properties of TB induced JJs under intrinsic chemical pressure V/x (created 
by the variation of the oxygen doping parameter 5). More specifically, a single SIS 
type junction (comprising a Josephson network) is formed around TB due to a lo- 
cal depression of the superconducting order parameter A(x) oc e(x) over distance do 
producing thus a weak link with (oxygen deficiency 5 dependent) Josephson coupling 
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J (5) = e(x)J = J Q (5)e-^l d ° where J (6) = e J = (fj, v /KQ )J (here J oc A /R n 
with R n being a resistance of the junction). Thus, the considered here model indeed 
describes chemically induced GS in underdoped systems (with 8 ^ 0) because, in accor- 
dance with the observations, for stoichiometric situation (when 8 ~ 0), the Josephson 
coupling J(8) ~ and the system loses its explicitly granular signature. 

To adequately describe chemomagnetic properties of an intrinsically granular su- 
perconductor, we employ a model of 2D overdamped Josephson junction array which 
is based on the well known Hamiltonian 

N N aa 

w = E^(!-««^)+E^ (!) 

ij ij ij 

and introduces a short-range interaction between Adjunctions (which are formed around 
oxygen-rich superconducting areas with phases <f>i(t)), arranged in a two-dimensional 
(2D) lattice with coordinates x ; = (xj,yj). The areas are separated by oxygen-poor in- 
sulating boundaries (created by TB strain fields e(Xy-)) producing a short-range Joseph- 
son coupling Jij = J (8)e~^ ij '' d . Thus, typically for granular superconductors, the 
Josephson energy of the array varies exponentially with the distance x^- = Xj — x 3 - 
between neighboring junctions (with d being an average junction size). As usual, the 
second term in the rhs of Eq.(l) accounts for Coulomb effects where qi = — 2erij is the 
junction charge with n.j being the pair number operator. Naturally, the same strain 
fields e(xy) will be responsible for dielectric properties of oxygen-depleted regions as 
well via the ^-dependent capacitance tensor Cij{8) = (7[e(xy)]. 

If, in addition to the chemical pressure V/x(x) = KQ Ve(x.), the network of super- 
conducting grains is under the influence of an applied frustrating magnetic field B, the 
total phase difference through the contact reads 

hj (t) = 0° + ^ (xy A n ?J ) • B + ^ n 13 , (2) 

where 0^- is the initial phase difference (see below), = Xy/| Xy | with = 
(xj + Xj)/2, and w = 2\ L (T) + I with X L being the London penetration depth of 
superconducting area and I an insulator thickness (which, within the discussed here 
scenario, is simply equal to the TB thickness [2l)]). 

To neglect the influence of the self-field effects in a real material, the corresponding 



Josephson penetration length Xj = J $o/2vr/xoj c w must be larger than the junction 
size d. Here j c is the critical current density of superconducting (hole-rich) area. As 
we shall see below, this condition is rather well satisfied for HTS single crystals. 

Within our scenario, the sheet magnetization M of 2D granular superconductor is 
defined via the average Josephson energy of the array 

r dt r d 2 x 

<H>= — Hhc,t) 3 

JO T J S 

as follows 

M(B,*) S -^>, (4, 
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where s = 2ird 2 is properly defined normalization area, r is a characteristic Josephson 
time, and we made a usual substitution ^ijif) — > - J d 2 xA(x,t) valid in the 

long- wavelength approximation \27\ . 

To capture the very essence of the superconducting analog of the chemomagnetic 
effect, in what follows we assume for simplicity that a stoichiometric sample (with 
8 ~ 0) does not possess any spontaneous magnetization at zero magnetic field (that 
is M(0, 0) = 0) and that its Meissner response to a small applied field B is purely 
diamagnetic (that is M(B,0) ~ — B). According to Eq.(4), this condition implies 

= 1-kvtl for the initial phase difference with m = 0, ±1, ±2, ... 

Taking the applied magnetic field along the c-axis (and normal to the CuO plane), 
that is B = (0, 0, B), we obtain finally 

M(B, 8) = -M (8) - ■ z ^ (5) 



'l + b 2 )(l + (b-b 



for the chemically-induced sheet magnetization of the 2D Josephson network. 

Here M (8) = J (8)/B with J (8) defined earlier (in what follows, M (0) is M (8 ~ 
0)), b = B/B , and b^ = BjB a ~ (k B Tr/h)8 where B^S) = (fi v T/h)B is the 
chemically-induced contribution (which disappears in optimally doped systems with 
8 ~ 0), and So = $o/wd is a characteristic Josephson field. Fig.l shows changes of the 
initial (stoichiometric) diamagnetic magnetization M/M (solid line) with oxygen defi- 
ciency 8. As is seen, even relatively small values of 8 parameter render a low field Meiss- 
ner phase strongly paramagnetic (dotted and dashed lines). The inset of Fig.l presents 
a true chemomagnetic effect with concentration (deficiency) induced magnetization 
M(0, 8) in zero magnetic field. According to Eq.(5), the initially diamagnetic Meissner 
effect turns paramagnetic as soon as the chemomagnetic contribution B^{8) exceeds an 
applied magnetic field B. To see whether this can actually happen in a real material, 
let us estimate a magnitude of the chemomagnetic field B M . Typically [13 ESI, f° r HTS 
single crystals Ax(0) ~ 150nm and d ~ lOnm, leading to -Bo — 0.5T. Using r ~ fij \x v 
and j c = 10 10 A/m 2 as a pertinent characteristic time and the typical value of the critical 
current density, respectively, we arrive at the following estimate of the chemomagnetic 
field B^S) ~ 0.5-Bq for 8 = 0.05. Thus, the predicted chemically induced PME should 
be observable for applied magnetic fields B ~ 0.5.Bo ~ 0.25T (which are actually much 
higher than the fields needed to observe the previously discussed [25 piezomagnetism 
and stress induced PME in high-T c ceramics). Notice that for the above set of param- 
eters, the Josephson length Aj ~ lfim, which means that the assumed in this paper 
small-junction approximation (with d Aj) is valid and the so-called "self-field" ef- 
fects can be safely neglected. So far, we neglected a possible field dependence of the 
chemical potential \i v of oxygen vacancies. However, in high enough applied magnetic 
fields B, the field-induced change of the chemical potential Afi v (B) = fi v (B) — (i v (0) 
becomes tangible and should be taken into account. As is well-known [2H1 121], m a 
superconducting state Afj, v (B) = —M(B)B/n, where M(B) is the corresponding mag- 
netization, and n is the relevant carriers number density. At the same time, within our 
scenario, the chemical potential of a single oxygen vacancy fi v depends on the concen- 
tration of oxygen vacancies (through deficiency parameter 8). As a result, two different 
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Figure 1: The magnetization M(B,5)/Mq(0) as a function of applied magnetic field 
B/Bq, according to Eq.(5), for different values of oxygen deficiency parameter: 5 ~ 
(solid line), 5 = 0.05 (dashed line), and 5 = 0.1 (dotted line). Inset: 5 induced 
magnetization M(0, 5)/M (0) in a zero applied magnetic field (chemomagnetism). 



effects are possible related respectively to magnetic field dependence of fi v (B) and to its 
dependence on magnetization fi v (M). The former is nothing else but a superconducting 
analog of the so-called magnetoconcentration effect (which was predicted and observed 
in inhomogeneously doped semiconductors [3111) with field-induced creation of oxygen 
vacancies c v (B) = c v (0) exp(—Afj, v (B) /ksT), while the latter (as we shall see in the 
next Section) results in a "fishtail" -like behavior of the magnetization. Let us start with 
the magnetoconcentration effect. Figure 2 depicts the predicted field- induced creation 
of oxygen vacancies 5(B) = 1 — c v (B) using the above-obtained magnetization M(B, 5) 
(see Fig.l and Eq.(5)). We also assumed, for simplicity, a complete stoichiometry of the 
system in a zero magnetic field (with 5(0) = 1 — c v (0) = 0). Notice that 5(B) exhibits 
a maximum at 5 C ~ 0.23 for applied fields B = Bq (in agreement with the classical per- 
colative behavior observed in non-stoichiometric Y Ba2Cu^O-j^s samples [T5 j l2*lj). 
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Figure 2: Magnetic field dependence of the oxygen deficiency parameter 5(B) (magne- 
toconcentration effect). 

Finally, let us show that in underdoped crystals the above-discussed osmotic mecha- 
nism of oxygen transport is indeed much more effective than a traditional diffusion. 
Using typical YBCO parameters [23], e = 0.01, Qo = a o with do = 0.2nm, and 
K = llbGPa, we have fi v (0) = €qKQq ~ lmeV for a zero-field value of the chemical 
potential in HTS crystals, which leads to creation of excess vacancies with concen- 
tration c v (0) = e -^(°)/ k BT 0.75 (equivalent to a deficiency value of 5(0) ~ 0.25) 
at T = T c , while the probability of oxygen diffusion in these materials (governed by 
a rather high activation energy Ud — leV) is extremely slow under the same condi- 
tions because D oc e~ Ud ' kBTc 1. On the other hand, the change of the chemical 
potential in applied magnetic field can reach as much as 29 j Afi v (B) ~ 0.5meV^ for 
B = 0.5T, which is quite comparable with the above-mentioned zero-field value of 
H v (0). Let us turn now to the second effect related to the magnetization dependence of 
the chemical potential [x v (M(B)). In this case, in view of Eq.(2), the phase difference 
will acquire an extra M(B) dependent contribution and as a result the r.h.s. of Eq.(5) 
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Figure 3: A "fishtail" -like behavior of magnetization m/ = M(B,5(B))/Mq(0) in 
applied magnetic field B/B in the presence of magnetoconcentration effect (with 
field-induced oxygen vacancies 5(B), see Fig. 2) for two values of field- free deficiency 
parameter: 5(0) ~ (solid line), and 5(0) = 0.1 (dashed line). 

will become a nonlinear functional of M(B). The numerical solution of this implicit 
equation for the resulting magnetization rrif = M(B, 5(B)) / Mq(0) is shown in Fig.3 for 
the two values of zero-field deficiency parameter 5(0). As is clearly seen, rrif exhibits 
a field-induced "fishtail" -like behavior typical for underdoped crystals with intragrain 
granularity (for symmetry and better visual effect we also plotted — rrif in the same 
figure). The extra extremum of the magnetization appears when the applied magnetic 
field B matches an intrinsic chemomagnetic field B^(5(B)) (which now also depends 
on B via the above-discussed magnetoconcentration effect). Notice that a "fishtail" 
structure of rrif manifests itself even at zero values of field-free deficiency parameter 
5(0) (solid line in Fig.3) thus confirming a field-induced nature of intrinsic granular- 
ity f3J EE 1221 123 El]- At the same time, even a rather small deviation from the 
zero-field stoichiometry (with 5(0) = 0.1) immediately brings about a paramagnetic 
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Meissner effect at low magnetic fields. Thus, the present model predicts appearance 
of two interrelated phenomena, Meissner paramagnetism at low fields and "fishtail" 
anomaly at high fields. It would be very interesting to verify these predictions exper- 
imentally in non-stoichiometric superconductors with pronounced networks of planar 
defects. 



3. MAGNETIC FIELD INDUCED POLARIZATION EFFECTS 

In this Section, within the same model of JJAs created by a regular 2D network 
of twin-boundary (TB) dislocations with strain fields acting as an insulating barrier 
between hole-rich domains in underdoped crystals, we discuss charge-related effects 
which are actually dual to the above-described phase-related chemomagnetic effects. 
More specifically, in what follows, we are interested in the behavior of magnetic field 
induced electric polarization (chemomagnetoelectricity) in chemically induced GS. 

Recall that a conventional (zero-field) pair polarization operator within the model 
under discussion reads [27J EH] 

N 

p = J2 ( 6 ) 

i=i 

In view of Eqs.(l), (2) and (6), and taking into account a usual "phase-number" commu- 
tation relation, rij] = i5ij, it can be shown that the evolution of the pair polarization 
operator is determined via the equation of motion 

= 4 [P; = T £ J « Sin ^'W^ ( 7 ) 
ij 

Resolving the above equation, we arrive at the following net value of the magnetic-field 
induced longitudinal (along a;-axis) electric polarization P(5, B) =< p x (t) > and the 
corresponding effective junction charge 

T t 2 

Q(S, B) = ^4 fdt fdt' f^- sin 0(x, f)xe-W\ (8) 
nrd J J J S 
o o 

where S = 2ird 2 is properly defined normalization area, r is a characteristic time (see 
Discussion), and we made a usual substitution j;J2ij ^-ijif) — ¥ § / d 2 xA(x., t) valid in 
the long- wavelength approximation |27j . 

To capture the very essence of the superconducting analog of the chemomagneto- 
electric effect, in what follows we assume for simplicity that a stoichiometric sample 
(with 5 ~ 0) does not possess any spontaneous polarization at zero magnetic field, that 
is P(0, 0) = 0. According to Eq.(8), this condition implies = lixm for the initial 
phase difference with m = 0, ±1, ±2, ... 

Taking the applied magnetic field along the c-axis (and normal to the CuO plane), 
that is B = (0, 0, B), we obtain finally 



9 



0.9 




Figure 4: The effective junction charge Q(8, B)/Q(8 C , 0) (chemomagnetoelectric effect) 
as a function of applied magnetic field B/Bq, according to Eq.(9), for different values 
of oxygen deficiency parameter: 5 ~ (solid line), 5 = 0.1 (dashed line), and 5 = 0.2 
(dotted line). 



for the magnetic field behavior of the effective junction charge in chemically induced 
granular superconductors. Here Qo(S) = erJ (5)/h with Jo(5) defined earlier, b = 
B/B , b = b — b^, and 6 M = B^/Bo ~ (kBTr/h)5 where B^S) = (fi v T/h)B is the 
chemically-induced contribution (which disappears in optimally doped systems with 
5 ~ 0), and B = § /wd is a characteristic Josephson field. 

Fig. 4 shows changes of the initial (stoichiometric) effective junction charge Q/Qo 
(solid line) with oxygen deficiency 8. Notice a sign change of Q/Qo (dotted and dashed 
lines) driven by non-zero values of S at low magnetic fields (a charge analog of chemically 
induced PME). According to Eq.(9), the effective charge changes its sign as soon as the 
chemomagnetic contribution B^S) exceeds an applied magnetic field B. At the same 
time, Fig.5 presents a true chemoelectric effect with concentration (deficiency) induced 
effective junction charge Q(S, 0) in zero magnetic field. Notice that Q(S, 0) exhibits 
a maximum around 5 C ~ 0.2 (in agreement with the classical percolative behavior 
observed in non-stoichiometric Y Ba^Cu^O-js samples |lfij). 

It is of interest also to consider the magnetic field behavior of the concomitant 
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Figure 5: Chemically induced effective junction charge Q(8, 0) / Q(S C , 0) in a zero applied 
magnetic field (true chemoelectric effect). 



effective flux capacitance C = rdQ(5, B)/d& which in view of Eq.(9) reads 



C(S,B) = C (S) 



1 - 366 - 36 2 + 66 3 
(1 + 6 2 )(1 + 6 2 ) 3 



(10) 



where <£> = SB, and C (5) = tQ (5)/$ . 

Fig. 6 depicts the behavior of the effective flux capacitance C(5, B)/Co(0) in applied 
magnetic field for different values of oxygen deficiency parameter: 5 ~ (solid line), 
5 = 0.1 (dashed line), and 5 = 0.2 (dotted line). Notice a decrease of magnetocapaci- 
tance amplitude and its peak shifting with increase of 5 and sign change at low magnetic 
fields which is another manifestation of the charge analog of chemically induced PME 
(Cf. Fig. 4). Up to now, we neglected a possible field dependence of the chemical po- 
tential fi v of oxygen vacancies. Recall, however, that in high enough applied magnetic 
fields B, the field-induced change of the chemical potential Afi v (B) = fi v (B) — /i v (0) 
becomes tangible and should be taken into account [2E1 E2] • As a result, we end up 
with a superconducting analog of the so-called magnetoconcentration effect with 
field induced creation of oxygen vacancies c v (B) = c v (0) exp(— A/x^(S)//c B T) which in 
turn brings about a "fishtail" -like behavior of the high-field chemo magnetization (see 
Section 2 for more details). Fig. 7 shows the field behavior of the effective junction 



11 




Figure 6: The effective flux capacitance C(5, B)/C(S C , 0) as a function of applied mag- 
netic field B/Bq, according to Eq.(10), for different values of oxygen deficiency param- 
eter: 5 ~ (solid line), 5 = 0.1 (dashed line), and 5 = 0.2 (dotted line). 

charge in the presence of the above-mentioned magnetoconcentration effect. As it is 
clearly seen, Q(5(B), B) exhibits a "fishtail" -like anomaly typical for previously dis- 
cussed [22] chemomagnetization in underdoped crystals with intragrain granularity (for 
symmetry and better visual effect we also plotted —Q(5(B),B) in the same figure). 
This more complex structure of the effective charge appears when the applied magnetic 
field B matches an intrinsic chemomagnetic field B tl (5(B)) (which now also depends on 
B via the magnetoconcentration effect). Notice that a "fishtail" structure of Q(5(B), B) 
manifests itself even at zero values of field-free deficiency parameter 5(0) (solid line in 
Fig. 7) thus confirming a field- induced nature of intrinsic granularity. Likewise, Fig. 8 
depicts the evolution of the effective flux capacitance C(5(B), B)/Cq(0) in applied mag- 
netic field B/Bq in the presence of magnetoconcentration effect (Cf. Fig. 6). Thus, the 
present model predicts appearance of two interrelated phenomena dual to the previ- 
ously discussed behavior of chemomagnetizm (see Section 2), namely a charge analog 
of Meissner paramagnetism at low fields and a charge analog of "fishtail" anomaly at 
high fields. To see whether these effects can be actually observed in a real material, let 
us estimate an order of magnitude of the main model parameters. 
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Figure 7: A "fishtail" -like behavior of effective charge Q(8(B), B)/Q(8 C , 0) in applied 
magnetic field B/B in the presence of magnetoconcentration effect (with field-induced 
oxygen vacancies 8(B)) for three values of field-free deficiency parameter: 8(0) ~ 
(solid line), 8(0) = 0.1 (dashed line), and 8(0) = 0.2 (dotted line). 

Using typical for HTS single crystals values of \l(0) ^ 150nm, d ~ 10nm, and 
j c ~ 10 10 A/m 2 , we arrive at the following estimates of the characteristic Bq ~ 0.5T 
and chemomagnetic B^(8) ~ 0.5i?o fields, respectively. So, the predicted charge analog 
of PME should be observable for applied magnetic fields B < 0.25T. Notice that, for 
the above set of parameters, the Josephson length is of the order of Aj ~ 1/im, which 
means that the assumed in this paper small-junction approximation is valid and the 
" self- field" effects can be safely neglected. 

Furthermore, the characteristic frequencies u ~ r _1 needed to probe the suggested 
here effects are related to the processes governed by tunneling relaxation times r ~ 
Ti/Jq(8). Since for oxygen deficiency parameter 8 = 0.1 the chemically-induced zero- 
temperature Josephson energy in non-stoichiometric YBCO single crystals is of the 
order of Jo(8) — ksTcS ~ lmeV, we arrive at the required frequencies of u ~ 10 13 Hz 
and at the following estimates of the effective junction charge Qo ~ e = 1.6 x 10~ 19 C 
and flux capacitance Cq ~ 10~ 18 F. Notice that the above estimates fall into the range of 
parameters used in typical experiments for studying the single-electron tunneling effects 
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Figure 8: The behavior of the effective flux capacitance C(5(B), B)/C(S C , 0) in applied 
magnetic field B/B in the presence of magnetoconcentration effect for three values 
of field-free deficiency parameter: 5(0) ~ (solid line), 5(0) = 0.1 (dashed line), and 
5(0) = 0.2 (dotted line). 



both in JJs and JJAs jT| 12 E3 EH] suggesting thus quite an optimistic possibility to 
observe the above-predicted field induced effects experimentally in non-stoichiometric 
superconductors with pronounced networks of planar defects or in artificially prepared 
JJAs. (It is worth mentioning that a somewhat similar behavior of the magnetic 
field induced charge and related flux capacitance has been observed in 2D electron 
systems 

And finally, it can be easily verified that, in view of Eqs.(6)-(8), the field-induced 
Coulomb energy of the oxygen-depleted region within our model is given by 

with Q(5,B) and C(5,B) defined by Eqs. (9) and (10), respectively. 

A thorough analysis of the above expression reveals that in the PME state (when 
B <C Bfj) the chemically-induced granular superconductor is in the so-called Coulomb 
blockade regime (with Ec > Jo), while in the "fishtail" state (for B > By) the en- 
ergy balance tips in favor of tunneling (with Ec < Jo). In particular, we obtain that 
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E C (S,B = O.lBp) = f J (S) and E C (6,B = B^) = § J (S). It would be also inter- 
esting to check this phenomenon of field-induced weakening of the Coulomb blockade 
experimentally. 



4. CONCLUSION 

In conclusion, within a realistic model of 2D Josephson junction arrays (created 
by 2D network of twin boundary dislocations with strain fields acting as an insulating 
barrier between hole-rich domains in underdoped crystals), a few novel effects expected 
to occur in intrinsically granular material are predicted, including phase-related and 
(dual) charge-related phenomena. The conditions under which these effects can be 
experimentally measured in non-stoichiometric high-T c superconductors were discussed. 
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